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Abstract: We compute the masses of the n and of the p mesons in the quenched 
approximation on a lattice with fixed lattice spacing a ~ 0.145 fm for SU(iV) gauge 
theory with N = 2,3, 4, 6. We find that a simple linear expression in 1/iV 2 correctly 
captures the features of the lowest-lying meson states at those values of N. This 
enables us to extrapolate to iV = oo the behaviour of function of the quark 

mass and of Trip £ls £l function of m n . Our results for the latter agree within 5% with 
recent predictions obtained in the AdS/CFT framework. 
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1. Introduction 

The theory of strong interactions, Quantum Chromodynamics (QCD), is a SU(3) 
gauge theory with n/ flavors of fermionic matter fields in the fundamental represen- 
tation of the color group. For sufficiently small n/, QCD displays many interesting 
non-perturbative phenomena, which are not captured by the conventional expansion 
in powers of the coupling constant. However, if we consider a SU(iV) gauge theory, 
with a generic number of colors N, in the limit where N becomes large, the pertur- 
bative expansion can be reorganized in powers of 1/N, and the contribution of each 
diagram can be directly related to its topology [p], ||. The leading contribution in 
this expansion is given by planar diagrams, and a simple power counting argument 
suggests that corrections are 0(1/ N 2 ) in a pure gauge theory, while the fermionic 
determinant yields corrections 0(nf/N). 

The large-iV expansion is a powerful tool to explore the strongly interacting 
regime of gauge theories, and recent developments in string theory have provided 
beautiful insights in our understanding of the planar limit through the gauge-gravity 
correspondence || (see 0] for an introductory review of recent developments). The 
lattice formulation of gauge theories allows one to study the non-perturbative dy- 
namics from first principles by numerical simulations, and can therefore be used to 
investigate how the N = oo limit is approached. A number of studies in recent 
years |6], [7|, |8|, ^, |I0J [TT] have analyzed in detail several features of pure gauge 



theories for N > 2, including the spectrum of glueballs, the fc-string tension, and 
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topology, both at zero and finite temperature. A very precocious scaling has been 
observed for all observables that have been considered so far, with 1/N 2 corrections 
being able to accommodate the values of the observables already for N = 3 and in 
most of the cases also for N = 2. 

The convergence to the large-iV limit for theories with fermions could also be 
addressed by dynamical simulations. The contributions of the fermionic determinant 
should increase the size of the corrections, as pointed out above. An intermediate 
step at a lesser computational cost is the study of properties of mesons and baryons 
in theories with quenched fermions. Note that since the fermionic determinant is 
suppressed in the large- N limit, simulations in the quenched approximation should 
converge to the same limit as in the theory with dynamical fermions, but with cor- 
rections 0(1/ N 2 ). 

This paper focuses on the low-lying states of the mesonic spectrum for SU(iV) 
theories in the quenched approximation and N = 2,3, 4, 6. By generalizing the lattice 
Dirac operator to handle spinors of arbitrary dimension in color space, we compute 
two-point functions for Wilson fermions at one value of the lattice spacing and several 
values of the bare quark mass. The mass dependence of the spectrum is studied, and 
extrapolated to the large-A limit. Our results are consistent with a 1/N 2 scaling, 
and the results for N = oo can be used as an input for analytical approaches that 
study the meson spectrum of strongly-interacting gauge theories. Some aspects of 
the meson spectrum at large N (in particular, the dependence of the mass of the 
pion from the quark mass) have also been investigated in fT2|j . 

We stress that this calculation is meant to be exploratory, trying to favor a 
first overall physical picture over more formal and technical points. A more detailed 
calculation is currently in progress and will be reported elsewhere. 

The paper is organized as follows. Sect, ^recalls the basic framework that is used 
for extracting the mesonic spectrum from field correlators in quenched lattice gauge 
theories and summarizes the choice of bare parameters for each value of N. The 
numerical results and their analysis are presented in Sect. |3]. Finally, we conclude 
by discussing the large- N extrapolation in Sect. f| and its relevance for AdS/QCD 
studies (see |13[ for a review) in Sect. ||. 



As this work was being completed we noticed that similar problems have been 
investigated in Ref. ||14j| . The preliminary results presented there are obtained on 
slightly smaller lattices (in physical units) with a finer lattice spacing. The two 
sets of data are complementary and in qualitative agreement. Future calculations 
will hopefully achieve precise continuum results for the large- N limit of the mesonic 
spectrum. 

2. Lattice formulation 

A Monte Carlo ensemble of gauge fields is generated using the Wilson formulation 



- 2 - 



of pure SU(iV) gauge theory on the lattice, defined by the plaquette action 

S = - A Tr [ U ( x > + ^ v)TP{x + p, /i)f/ f (x, v) + h.c] , (2.1) 

X,)X>V 

where U(x,fi) G SU(iV) are the link variables. The link variables are updated using 



a Cabibbo-Marinari algorithm fT5| , where each SU(2) subgroup of SU(iV) is updated 
in turn. We have alternated microcanonical and heat-bath steps in a ratio 4:1. We 
call sweep the sequence of four microcanonical and one heat-bath update. 
The action of the massive Dirac operator on a generic spinor field 4>(x) is: 

D m (j){x) = (D + m)<j)(x) 

= -- I ^2 [(1 - 7 M ) U (x, n)(j){x + fj) + (1 + 7^) U(x - //, fiY4>(x - n ) 

-(8 + 2m)0(x)}. (2.2) 

The bare mass is related to the hopping parameter used in the actual simulations by 

l/(2«) = 4 + m. (2.3) 

The complete set of bare parameters used in our simulations is summarized in Tab. [I]. 
For this preliminary study we have run simulations for values of N ranging from 2 
to 6 at one value of the lattice spacing only. The values of (3 have been chosen in 
such a way that the lattice spacing is constant across the various N. More in detail, 
at each N we chose the critical value of (3 for the deconfinement phase transition 
at ^ = 5 [0]. To set the scale, another physical quantity (like e.g. the string 
tension) can be used. Since different quantities have different large-iV corrections, 
a different choice for the scale will affect the size of the 1/N 2 corrections, but not 
the N = oo value. Using the value T c = 270 MeV, for the lattice spacing a we get 
a ~ 0.145 fm. For our simulations, we used a x N t lattice with the spatial size 
N s = 16 (which corresponds to about 2.3 fm) and the temporal size N t = 32 (about 
4.6 fm in physical units). For the quark masses used in this work, our calculation 
should be free from noticeable finite size effects. For the fermion fields we used 
periodic boundary conditions in the spatial directions and antiperiodic boundary 
conditions in the temporal direction. For the gauge field we used periodic boundary 
conditions in all directions. 

We have performed a chiral extrapolation using data from meson correlators at 
five values of k for each N. The choice of values for k relies on previous experience 
with SU(3) simulations to yield pseudoscalar meson masses > 450 MeV. The same 
values of k have been used for all values of iV > 3, while for SU(2) a different choice 
turned out to be necessary, since all k's but the lowest one were higher than k c . 
Because of the different additive renormalization, these values of k yield different 
values for the bare PCAC mass as iV is varied. 
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N 


P 


K 


2 


2.3715 


0.156, 0.155, 0.154, 0.153, 0.152 


3 


5.8000 


0.161, 0.160, 0.159, 0.1575, 0.156 


4 


10.6370 


0.161, 0.160, 0.159, 0.1575, 0.156 


6 


24.5140 


0.161, 0.160, 0.159, 0.1575, 0.156 



Table 1: Bare parameters used in our simulations. 



The mesonic spectrum is extracted from the zero-momentum two-point correla- 
tors of quark bilinears with the quantum numbers required to interpolate between 
a meson state and the vacuum. Let I\ and T 2 be two generic products of Dirac 7 
matrices, a two-point correlator is defined as 

CWt) = Yl ((^1^ (*> x ) (uT 2 d) (0)} , (2.4) 

X 

where u and d are the fields corresponding to two different quark flavors, which 
from now on we take to be mass degenerate. All possible choices for the non-singlet 
quark bilinears and the quantum numbers of the corresponding physical states are 
summarized in Tab. 



Particle 


a 


7T 


p 




h 


Bilinear 


ud 


u^d, U'jo^d 






uji^/jd 


jPC 


0++ 


0-+ 


\ — 


1++ 


1+- 



Table 2: Bilinear operators for the computation of non-singlet meson masses. 



CWt) = ~J2\ tT [7or t 1 7oG(x,0)r 275 G'(x,0) t 75 

X 

= - ( tr UoThoH(x, 0h 5 T 2 H(x, 0) 75 



Performing the Wick contractions we can rewrite Cr 1 ,r 2 (^) m terms of the quark 
propagator G(x, y) = (D m ) _1 (s, y) or, equivalently, of its hermitean version H(x, y) = 
G(x,y)^ 5 : 



(2.5) 



(tr indicates the trace over spinor and color indices). 

The propagator G(x, 0)ab is obtained by inverting the Dirac operator D m over 
point sources (capital roman letters A, B, . . . are used for collective indices over spin 
and color): 



G{x, 0)ab = (D m ) A c(x,y)8 C BSy,o = {D m ) A l c {x,y)rfc\y) 
where the second equality defines the AN point sources rj( B \ 



(2.6) 
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The algorithm used for the inversion in Eq. Q2.6Q , is a multishift QMR. This 
enables us to compute all the quark propagators corresponding to different masses 
simultaneously. We use a version of the QMR suitable for 75-hermitean matrices 
with even-odd preconditioning of the Dirac matrix In the rare cases when the 



algorithm fails to converge, we continue the search for a solution using the MINRES 
algorithm, which is guaranteed to converge, on the hermitean version of the Dirac 
operator. For all the inversions we required a relative precision of 10 -5 . To this 
accuracy, and for the parameters given above, the number of required applications 
of the Dirac operator to compute the propagator G(x, 0)ab at fixed values for the 
hopping parameters k is found to become independent of N. We found that the 
average number of applications of the Dirac matrix required is about 7500, 5000, 
5000 for iV = 3,4,6 respectively (for SU(2) we used different parameters). 

From general large-iV arguments, we expect the occurrence of exceptional con- 
figurations to be suppressed as N increases. This is confirmed by preliminary results 
reported in Ref . ]TJ| . At the values of the parameters we have simulated, there is no 



sign of the presence of exceptional configurations. 

Simulations have been performed with a bespoke code, which has been tested 



against published results for SU(3) (see |17| for a review of the literature). We have 



collected 100 configurations for each value of N, separated by 50 Monte Carlo sweeps. 

3. Numerical results 

3.1 Extracting masses from correlators 

Masses can be extracted from the large-t behavior of Cr,r{t)- Inserting the energy 



eigenstates in the RHS of Eq. ( |2.4| ) yields 

Cr,r(t) = $> 4 | 2 e-"^ , (3.1) 

i 

where c$ = (l/2m;) (0 \ (uTd) (0)| i), T is one of the 7 matrix products appearing in 
the bilinears in Tab. 0, \i) is an eigenstate of the Hamiltonian with the same quantum 
numbers as the fermion bilinear, and is the mass of the \i) eigenstate. In the limit 
t — > 00, the previous equation becomes 

Cr,r(t) = \co\ 2 e- mot , (3.2) 

t— >oo 

i.e. at large time correlation functions decay in time as a single exponential with 
a typical time given by the inverse mass of the lowest-lying state in the spectrum 
with matching quantum numbers. The lowest mass in a given channel can then be 
extracted as 

m = — lim ■ . (3.3) 

t^oo t 
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t 



Figure 1: Correlators for SU(4) and T = 75 at the values of k shown. The correlators 
have been normalized in such a way that at t = their value is 1. The lines joining the 
data are only guides for the eyes. 

Practically, one defines the effective mass mo(t) as 

m M = -**<g0T) (3 ' 4) 

and obtains m by fitting mo(t) to a constant at large enough t. 

On a finite lattice the exponential in the large-time behavior of the propagator 
is replaced by a cosh and an effective mass can be defined as 



(3.5) 



Typical examples of correlation functions and effective masses as a function of the 
separation between source and sink are shown respectively in Fig. p] and Fig. ||. 

We have estimated the errors on the correlators using a jack-knife method, and 
checked that the bootstrap method gives similar results. With simple link operators, 
we have been able to extract an unambiguous signal for the n and the p mesons (to 
which we limit our analysis). Other correlators yielded a noisy signal, and we plan to 
investigate the possibility of improving the signal-to-noise ratio by more sophisticated 
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Figure 2: Effective masses from the correlators in Fig. [T|. The straight lines are fits to the 
data at plateau. 



measurements. Due to the small number of data points, often correlated fits proved to 
be unreliable, as already observed (see e.g. Refs. |Tj| |19|). Hence, masses have been 
extracted with uncorrelated fits and the error estimated with a jack-knife procedure. 
We have checked that the uncorrelated fit results coincide with the correlated fit 
results whenever the correlated fits give reasonable values for the parameters and 
the x 2 - 

Our results for the various N for the PCAC mass mpcAc? the mass of the pion 
m n , and the mass of the p m p are reported in Tabs. 0-^. The details of our analysis 
are explained in the following two subsections. 

In order to convert the results expressed in lattice units to masses in physical 
units, we note that 

am = (aT c )(m/T c ) = m/(5T c ) (3.6) 

(the last equality uses the fact that the lattice spacing has been fixed in such a way 
that the deconfinement phase transition corresponds to N t = 5). As a reference scale, 
we can use T c for SU(3), which is approximately 270 MeV. 

3.2 Meson masses at finite N 

The pion is the would-be Goldstone boson of chiral symmetry breaking. Chiral 
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K 


amp C Ac 


^pcac (MeV) 


am^ 


m u (MeV) 


am n 


m p (MeV) 


0.152 


0.0824(30) 


111(4) 


0.541(3) 


730(4) 


0.620(6) 


837(8) 


0.153 


0.0669(28) 


90(4) 


0.492(4) 


664(5) 


0.584(7) 


788(9) 


0.154 


0.0522(25) 


70(3) 


0.441(4) 


595(5) 


0.547(7) 


738(9) 


0.155 


0.0396(24) 


53(3) 


0.389(4) 


525(5) 


0.510(10) 


789(13) 


0.156 


0.0261(22) 


35(3) 


0.319(7) 


430(9) 


0.458(17) 


618(23) 



Table 3: Numerical results for SU(2). Masses in lattice units have been converted to 
physical units by noting a = 1/(5T C ). 



K 


ampcAc 


mpcAc (MeV) 


am^ 


(MeV) 


arrip 


m (MeV) 


0.156 


0.1047(24) 


141(3) 


0.625(2) 


844(3) 


0.720(3) 


972(4) 


0.1575 


0.0797(21) 


108(3) 


0.553(2) 


747(3) 


0.667(4) 


900(5) 


0.159 


0.0574(17) 


77(2) 


0.476(2) 


643(3) 


0.616(5) 


832(7) 


0.160 


0.0431(16) 


58(2) 


0.420(2) 


567(3) 


0.582(6) 


786(8) 


0.161 


0.0299(14) 


40(2) 


0.362(3) 


489(4) 


0.550(7) 


743(9) 


Table 4: 


Numerical results for SU(3). 


Masses in 


lattice units have been 


converted to 


physical units by notin 


go = l/(5T c ). 










K 


ampcAC 


mpcAc (MeV) 


am^ 


m w (MeV) 


am p 


m p (MeV) 


0.156 


0.1506(35) 


203(5) 


0.733(1) 


990(1) 


0.817(2) 


1103(3) 


0.1575 


0.1234(29) 


167(4) 


0.667(1) 


900(1) 


0.766(2) 


1034(3) 


0.159 


0.0981(22) 


132(3) 


0.598(1) 


807(1) 


0.714(2) 


964(3) 


0.160 


0.0817(19) 


110(3) 


0.549(2) 


741(3) 


0.680(2) 


918(3) 


0.161 


0.0659(17) 


89(2) 


0.499(2) 


674(3) 


0.646(3) 


872(4) 


Table 5: 


Numerical results for SU(4). 


Masses in 


lattice units have been 


converted to 



physical units by noting a = 1/(5T C ). 



perturbation theory at leading order predicts 



I 1 x ^ 



m n = A I . (3.7) 

\k k c J 

Hence the value of k corresponding to the chiral limit, be obtained by fitting 

the pion mass according to Eq. (|3.7|) . Eq. ( p.7|) is modified for the quenched theory, 
where quenched chiral logs appear. For quenched SU(iV) gauge theory we expect 

/j 1 \ V[2(l+«5)] 

m n = A( , (3.8) 

where 5 is positive, O(10 _1 ) for SU(3) and goes like 1/N [ BD] . 
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Figure 3: Effective masses from C 75i75 at k = 0.156. 



K 


ampcAC 


rapcAc (MeV) 


am-x 


m n (MeV) 


arrip 


m p (MeV) 


0.156 


0.1789(22) 


241(3) 


0.814(1) 


1099(1) 


0.889(2) 


1214(3) 


0.1575 


0.1509(20) 


203(3) 


0.752(1) 


1015(1) 


0.838(2) 


1131(3) 


0.159 


0.1248(18) 


168(2) 


0.687(1) 


927(1) 


0.786(2) 


1061(3) 


0.160 


0.1078(17) 


146(2) 


0.6425(9) 


867(1) 


0.752(2) 


1015(3) 


0.161 


0.0914(15) 


123(2) 


0.5952(9) 


804(1) 


0.718(2) 


969(3) 



Table 6: Numerical results for SU(6). Masses in lattice units have been converted to 
physical units by noting a = 1/(5T C ). 



The mass of the pion can be extracted by looking at correlators Cr,r i n which F 
is either 75 or 7075. In the latter case, it was not possible to extract a signal for all 
k's in SU(2). Hence, although in general mass plateau fits of 7075 have a lower % 2 , 
we will mostly focus on numerical results obtained with T = 75. 

As N grows, so does m n at fixed k. A plot comparing numerical results for 
the effective mass extracted from C 75i75 is shown in Fig. |3|. A linear fit to the data 
according to Eq. ( |3.7| ) enables us to extract the critical values of k for the at 
which we have simulated. Results are shown in Tab. [7[ The mass obtained from the 
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N 




A 


2 


0.15827(12) 


1.0583(99) 


3 


0.16359(28) 


1.142(21) 


4 


0.16556(23) 


1.201(14) 


6 


0.16716(12) 


1.2422(69) 



Table 7: Fitted values for k c and A at various N. 



C 7o75i7075 correlator yields compatible results. Higher statistics and a careful study 
of the systematics are necessary for a more precise determination of the critical value 
of K. 

Generally the reduced x 2 of the fits according to Eq. (|3.7| ) varies between two and 
four. One can check whether this relatively high value of \ 2 r is due to the fact that we 
are neglecting chiral logarithms. For N >2, fits according to (p.8|) yield a value of xl 



that is below one, but S is found to be negative. This agrees with the findings of [|T7 
where the negative value is interpreted as a consequence of simulating far from the 
chiral limit. In fact, for N = 2, where we have the lightest pion mass, 5 is found to 
be compatible with zero. As one would have expected, for m n > 450 MeV there is no 
sensitivity to the chiral logarithms [0. Instead of relying on phenomenological fits 
like in [[17]], we acknowledge the impossibility to determine 5 and neglect the chiral 
logarithms, using the chiral behavior ( |3.8| ) to get an estimate for the systematic error 
associated with this approximation. For SU(2), the three-parameter fit gives a value 
of k c that is higher than the fit with 5 = by about 1%. Considering that 5 oc N" 1 , 
a conservative but safe estimate for the systematic error associated with the chiral 
logarithms is of the order of a few percent. This is in agreement with the literature 
for SU(3) [17]]. We will come back on issues associated with the chiral logarithms in 
the next subsection. 

The mass of the p has been extracted from C 7ij7i , after taking the average over 
the spatial direction i of the correlation functions. Fits in the C 707ii707i channel also 
yield compatible results. 

At small quark mass, m p depends linearly on the quark mass and goes to a finite 
value in the chiral limit. Using Eq. ( |3.7| ), this can be rephrased into the following 
relationship between m p and m n 

m p = m x p + Bml , (3.9) 

where is the mass of the p meson at the chiral point. Note that the previous 



equation is not modified by chiral logarithms [20|| . Assuming that Eq. ( |3.9| ) holds in 
our case 1 , we can fit and B at the various values of N from our data. Our results 



1 More sophisticated dependencies (e.g. the addition of a linear term in to Eq. (3J5), which 
is motivated by phenomenology) are also supported by our data. In the absence of any evidence 
against it, we chose to fit the parameters using the simple chiral functional behavior. 
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N 


am* 


B 


2 


0.3890(75) 


0.797(31) 


3 


0.4683(25) 


0.6455(84) 


4 


0.5018(36) 


0.5905(88) 


6 


0.5238(40) 


0.5533(77) 



Table 8: Extrapolation of m p to the chiral limit. 



for those quantities are reported in Tab. |8]. The reduced \ 2 °f the fits (which keep 
into account both the error on m p and the error on m n ) is always less than one. 

3.3 PCAC 

As noted in the previous subsection, Eq. ( |3.7| ) only holds for the full theory, while 
it is modified at small masses, where quenched chiral logs become important. An 
alternative way of defining k c , which is free from these ambiguities, makes use of 
the partially conserved axial current (PCAC) relation. In the continuum, the PCAC 
relation reads 

d^(x) = 2m PCAcJ (x) , (3.10) 

with A^(x) = u{x)^/y^d{x) and j = u(x) r y 5 d(x) . The previous equation allows us to 
determine mpcAc as 

l (Jdf(3„/l°(x))jt(y))) 

mpcAC " 2 <Jdfi(*)j%» ' (3 U) 

where y is an arbitrary point. On the lattice an effective mass mpc\c{t) can be 
defined as 

A/A ^ ^7075 ,75 (j ~j~ jj ~ ^7075,75 ~ jj /o in\ 

and once again fitted at plateau. Note that with our choice for the discretized 
fermions PCAC holds on the lattice up to terms 0(a). In practice, since ^pcac(^) 
defined through Eq. Q3.12j ) is antisymmetric around the point N t /2, one averages 
the absolute values at points t and N t — t. An example of an effective mass plateau 
obtained using Eq. ( |3.12| ) is given in Fig. 

Since m PC Ac = Z m (l/n — l/« c ), we can determine k c as the value for which 
^pcac = 0. A linear fit to the data enables us to extract k c . Results for = 2, 3, 4, 6 
are reported in Tab. f| Comparing with the similar fits from m n (Tab. 0), it is 
immediate to see that using m PC Ac we get values for k c that are systematically 
lower 2 . Although this effect is below half a percent, it is by far larger than the 

2 This should be contrasted with fits that keep into account chiral logarithms, for which we find 
values of k c systematically higher. 
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Figure 4: mpcAc(i) as a function of t for k = 0.156. 



N 




A 


2 


0.15792(25) 


0.331(15) 


3 


0.16306(10) 


0.372(10) 


4 


0.16513(15) 


0.422(12) 


6 


0.16657(15) 


0.438(11) 



Table 9: Extrapolation of n to the chiral limit using mpcAC- 



statistical errors. This discrepancy might be due to the different chiral behavior of 
the two definitions of the quark mass for the quenched theory or be a consequence 
of the underestimation of the errors due to the use of uncorrelated fits, as we have 
discussed above. However, at our value of the lattice spacing discretization errors 
also play a relevant part. In order to investigate these issues, we have analyzed 
as a function of m PC Ao Our results are reported in Fig. |5], where m n is plotted as a 
function of mpcAC- 

The expected quadratic behavior 

ml = CVipcac (3.13) 
is not obeyed by our data. To extrapolate to the chiral limit, we need to correct the 
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Figure 5: m 2 as a function of topcac at the various N. The values of m, at mpcAC = 
have been obtained with a linear fit to the data, as discussed in the text. 

above relationship by allowing for a non-zero value for m n when mpcAC is zero: 

ml = Cmp CAC + B , (3.14) 

where B and C depend on N. B (obtained from a fit according to the previous 
equation keeping into account both the errors on and on mpcAc) is roughly of 
order 10~ 2 and independent of N. If the existence of a constant (as a function of N) 
residual m n asm->0 were a sign of the failure of the quenched approximation, we 
would have expected B to go to zero as N — > oo. This expectation is not supported 
by our data. On the other hand, having fixed the lattice spacing across the gauge 
groups, any discretization artifact would be constant in N. Hence, it is likely that 
this residual mass is (mostly) due to the violation of PCAC on the lattice. If this is 
correct, also the systematic discrepancy between the two sets of n c should be due to 
lattice artifacts. In order to settle this issue, a systematic study at different lattice 
spacings needs to be performed. 



4. Extrapolation to SU (oo) 

Using data at finite N, we can estimate the behavior of the lowest-lying meson masses 
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at iV = oo. Following similar analysis performed in pure gauge ||, §, [7], ||, |], 1(| |TT|], 
we use predictions from the large- N expansion to see whether they hold in the non- 
perturbative regime. In practice, we take the asymptotic expansion for an observable 
O in the quenched case JTJ 

0(N) = 0(oo) + J2^ (4-1) 

i 

and we check whether a reasonable (as dictated by the number of data) truncation 
of this series accommodates our numerical values. In the pure gauge case, a preco- 
cious onset of the large-iV behavior has been found for all the observables that have 
been studied (which include glueball masses, deconfming temperature and topolog- 
ical susceptibility): the 0(1/N 2 ) correction correctly describes the data down to at 
least N = 3, often including also the case N = 2. From a qualitative point of view, 
it is already clear from what we have seen so far that the quantities we have investi- 
gated have a mild dependence on N. In this section, we want to study whether this 
dependence is correctly described by a large- iV-inspired expansion. 

k c can be computed in lattice perturbation theory [PT]] . The result at one loop 
is in agreement with the predictions of the large-iV limit: this quantity receives a 
correction 0(1/N 2 ). This motivates the fit 

k c (N) = Kc (oo) + ^ . (4.2) 

For k c obtained via Eq. fl3.7|), we get k c (oo) = 0.1682(1) and a = —0.0398(6), with 
X 2 . = 0.6. The quality of the fit is good, and the coefficient of the 1/N 2 correction 
is small, as one would expect for a series expansion. Similarly to the pure gauge 
case, we observe an early onset of the asymptotic behavior, which captures also the 
SU(2) value. Our data and the large-iV extrapolation are plotted in Fig. || The 
same extrapolation for the critical value of k obtained using the PCAC relation 
yields k c (oo) = 0.1675(2), a = —0.039(1) and Xr — 1-3- The discrepancy between 
the values of kJoo) could be due to lattice discretization artifacts, as discussed in 



Sect. |373|. We take the difference between the two determinations should be seen as 
an estimate of the systematic error. The fact that the angular coefficient a has the 
same value seems to corroborate this hypothesis. A more precise determination of 
n c is beyond the scope of this work. 

The slope A in Eq. (|3.7|) can also be extrapolated to the N = oo limit, with 



corrections that are 0(1/ N 2 ): 

A{N) = A(oo) + a/N 2 . (4.3) 

For the C 75>75 results, the fit gives A(oo) = 1.262(8) and a = -0.82(6) with xl = 1-2. 
This allows us to write the mass of the 7r as a function of k at iV = oo as 

= 1.262(8) (1/k - 5.945(4)) 1/2 , (4.4) 
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Figure 6: Extrapolation of k c to N = oo. 
where the values of k c obtained from fits to C 75 75 have been used (using the PCAC 



value gives a slightly discrepant result, for the reasons discussed in Sect. |3.3|) . We plot 
in Fig. the data for the dependence of m% as a function of k, with a fit according 
to Eq. (|3.7|). We also plot in the same figure Eq. Q4.4|). 

The parameters describing function of m n (Eq. (|3.9|)) are also expected 

to follow the asymptotic expansion (|4.1| ). A fit with only the leading 1/N 2 correction 
yields 

m*(A0 = 0.539(3) - 0.62(3)/iV 2 , xl = 0.008 ; 

B(N) = 0.5224(8) + 1.10(l)/iV 2 , xl = 0.7 . 1 ' J 

The tiny for B(oo) is particularly surprising, given the statistical independence of 
the measured values of B(N). We see that once again the leading behavior describes 
very well the parameters and that the coefficient of the 1/N 2 correction is order one. 
As a result of this analysis, at iV = oo we can describe function of m, 

as 

m p = 0.539(3) + 0.5224(8)m 2 . (4.6) 
This relationship, together with the data and the fits at finite N, is plotted in Fig. Rl 
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Figure 7: m n as a function of 1/k. The curves through the data are obtained from a fit 
assuming the expected leading dependence from chiral perturbation theory. Also shown is 
the extrapolation to N = oo. 

5. Discussion and conclusions 

In this paper, we have used standard lattice QCD methods for computing correla- 
tion functions to extract the masses of the p and tt mesons in the large- N limit of 
SU(iV) gauge theories. The masses in the limiting theory have been obtained with an 
extrapolation from the quenched data at N = 2,3,4,6 using the large-iV behavior 
deduced from arguments inspired by a diagrammatic expansion. We find that the 
extrapolation works well in its simplest form, i.e. using only the leading correction to 
the large- N value. This allows us to determine the behavior at N = oo of the mass 
of the pion as a function of the renormalized quark mass and of the mass of the p as a 
function of the mass of the pion (chiral perturbation theory has also been used as an 
input). The two central results of this paper are summarized by the parametrization 

mAN) = (062,8, - g - 5.945,4, + * 

and 

m p (N) = f 0.539(3) - °-^) + ((15224(8) + ml(N) . (5.2) 
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Figure 8: m p vs. m w . Lines through the data have been obtained with a plot inspired by 
chiral perturbation theory. Also shown is the extrapolation to N = oo. 



One of the motivations to perform a calculation from first principles for the 
large- A" limit of SU(iV) gauge theories is to compare with predictions obtained in 
the AdS/CFT framework. To date, no AdS background has been found that can 
be considered a good dual description of non-supersymmetric QCD extrapolated 
at large-A^. Hence, one could use our calculation to benchmark the proposed AdS 
ansatz. From this point of view, we notice a striking agreement with a calculation [E2 



13fl using the Constable- Myers background ||23|| , which (after normalizing the mass 



of the p in the chiral limit to our data) finds for the coefficient of (see Eq. ( |3.9|) ) 
0.57. This number is in agreement with our calculations within 5%. However, before 
we can draw any conclusion from the comparison with the lattice, an extrapolation 
of the lattice data to the continuum limit is needed. 

Beyond the specific numerical details, our calculation seems to indicate that (a) 
the large- N theory is a well defined theory; (b) lattice calculations can be successfully 
exploited to compute the parameters of this theory; (c) at least in the quenched 
theory, to describe results at any finite A" only the first term in the expected power 
series in 1/A^ 2 is required; (d) the coefficient of the correction is at most order one, 
justifying the idea of a power expansion. All these indications are perfectly in line 
with what we have already learned for the SU(AQ theory without fermionic matter. 
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Although this can be considered obvious, since our calculation is quenched, we stress 
that the N = oo limit is also quenched. In other words, to describe the limiting 
theory a quenched calculation suffices. The inclusion of the full fermion determinant 
becomes mandatory if one is interested in the actual size of the finite N corrections. 
In particular, one expects larger corrections (0(rif/N)) in the unquenched theory. 

As we have stressed several times, one of the main limitations of this calculation 
is that our chiral extrapolations are not sensitive to the expected chiral log behavior. 
We have conservatively estimated that this approximation produces a 3% systematic 
error. This error does not affect our conclusions. Moreover, we note that chiral 
logarithms do not modify Eq. (|5.2| ). In any case, in order to make more robust 
estimates, better control on the chiral extrapolation should be achieved. This requires 
simulating at smaller pion masses. It would also be nice to check that the chiral log 
effects decrease as N increases. 

Another source of systematic error in our calculation is the fact that simulations 
have been performed at one single lattice spacing. For this reason we regard our 
results as exploratory. As discussed in Sect. |3~3| , lattice artifacts do seem to play a 
role, although they are not big enough to spoil the features of the theory and the 
way in which the large-iV behavior is approached. Nevertheless, a study closer to 
the continuum is necessary to clarify these issues. Work for the extrapolation to the 
continuum limit is already in progress. For this extrapolation, the use of an improved 
fermion action can mitigate the discretization artifacts, decreasing considerably the 
required numerical effort. We shall explore this possibility in the future. As for finite 
size effects, we have argued in Sect. [2] that with our choice of parameters they can 
be neglected, but this also ought to be verified directly. 

Aside from technicalities, other features of the large-iV theory, like the spec- 
trum in the flavor singlet channel and masses of heavier mesons, also deserve to be 
investigated. While the latter problem can probably be dealt with using improved 
techniques for computing correlation functions (e.g. with smeared links replacing 
straight links and smeared sources replacing point sources), for the scalar mesons, 
for which disconnected contributions are important, different strategies need to be 



adopted. An adaptation of the techniques exposed in |24], ^] is in progress. Results 
will be reported in a future publication. 
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